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Abstract 

Griffin, Wu and Stringari have derived the hydrodynamic equations of a 
trapped dilute Bose gas above the Bose-Einstein transition temperature. We 
give the extension which includes hydrodynamic damping, following the clas- 
sic work of Uehling and Uhlenbeck based on the Chapman-Enskog procedure. 
Our final result is a closed equation for the velocity fluctuations <5v which 
includes the hydrodynamic damping due to the shear viscosity rj and the ther- 
mal conductivity k. Following Kavoulakis, Pethick and Smith, we introduce 
a spatial cutoff in our linearized equations when the density is so low that 
the hydrodynamic description breaks down. Explicit expressions are given 
for r\ and k, which are position-dependent through dependence on the local 
fugacity when one includes the effect of quantum degeneracy of the trapped 
gas. We also discuss a trapped Bose-condensed gas, generalizing the work of 
Zaremba, Griffin and Nikuni to include hydrodynamic damping due to the 
(non-condensate) normal fluid. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db 
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I. INTRODUCTION 



Recent observation of Bose-Einstein condensation in trapped atomic gases has stimulated 
interest in the collective oscillations of non-uniform Bose gases The hydrodynamic 

equations of a trapped Bose gas (neglecting damping) have been recently derived above 
H and below M the Bose-Einstein transition temperature (Tbec)- These hydrodynamic 
equations describe low frequency phenomena (uj <C l/r c , where r c is the mean time between 
collisions of atoms), when collisions are sufficiently strong to ensure local thermodynamic 
equilibrium. These results can be used to derive closed equations for the velocity fluctuations 
in a trapped gas, whose solution gives the normal modes of oscillation. In particular, above 
Tbec, Ref. || gives explicit results for the frequencies of surface, monopole, and coupled 
monopole-quadrupole modes. Available experimental data on the frequencies and damping 
of collective oscillations in trapped gases is mainly for the collisionless regime, rather 
than the hydrodynamic regime we consider in this paper. However sufficiently high densities 
which allow one to probe the hydrodynamic region at finite temperatures are expected in 
the near future (see also Refs. []5]||). Thus it is worthwhile to extend the discussion of 
hydrodynamic modes in Refs. HQ to include damping. 

In the present paper, we do this by following the standard Chapman-Enskog procedure, as 
first generalized for quantum gases by Uehling and Uhlenbeck [|7|,[|. In Section II, we derive 
the hydrodynamic equations of a trapped Bose gas above T B ec which include damping due 
to shear viscosity 77 and thermal conductivity k. We also obtain integral equations which give 
expressions for the transport coefficients i] and k, and in Section III we solve for these. Our 
approximation corresponds to a lowest-order polynomial approximation of the Chapman- 
Enskog results such as obtained by Uehling ||. The transport coefficients have a position 
dependence when we include the effect of quantum degeneracy of the trapped Bose gas. 

In Section IV, we use our results to discuss the hydrodynamic damping of the surface 
normal modes in a trapped Bose gas discussed in Ref. ||. In the low density tail of a trapped 
gas, the hydrodynamic description ceases to be valid. It is crucial to include a cutoff in the 
linearized hydrodynamic equations to take this into account, as pointed out by Kavoulakis, 
Petthick and Smith MM. We also explicitly show that with this spatial cutoff, the damping 
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of normal modes given by our linearized hydrodynamic equations agrees with the expression 
on which Ref. is based. 

In Section V, we extend these results to the Bose-condensed region just below T B ec- 
Zaremba, Griffin and Nikuni 0] have recently given an explicit closed derivation of two-fluid 
hydrodynamic equations starting from a microscopic model of a trapped weakly-interacting 
Bose-condensed gas. We generalize the results of Ref. || to include hydrodynamic damping. 
Finally in Section VI, we give some concluding remarks. 

II. CHAPMAN-ENSKOG METHOD 

In Section II- IV, we limit ourselves to a non-condensed Bose gas above T BEC , as in 

Ref. [Q. The atoms are described by the semi-classical kinetic equation for the distribution 

function /(r, p, t) of a Bose gas |L0| 
d p 



+ £:• V r - VU(r,t) ■ V p 



dt m 



/(r,M=f 



(1) 



coll 

Here U(r,t) = Uo(r) + 2gn(r,t) includes an external potential Uo(r) as well as the Hartree- 
Fock (HF) self-consistent mean field 2gn(r,t), where n(r,t) is the local density. In usual 
discussions @,|], the HF field is omitted but it will play a crucial role when we discuss a 
Bose-condensed gas below T BEC (see Section V). The quantum collision integral in the right 
hand side of (P is given by ]10) 

dt 



coll 



Ang j (2^1 (^y (2^ 5(p+p2 - p3 - p4) H^ + ^ - ^~^J 

x[(l + /)(1 + f 2 )f z f A ~ // 2 (1 + + A)], (2) 

where / = /(r, p, t),fi = /(r, Pi,t). The interaction strength g = 4irh 2 a/m is determined 
by the s-wave scattering length a. 

The conservation laws are obtained by multiplying (1) by 1, p and p 2 and integrating 
over p. In all three cases, the integrals of the collision term in (H) vanish and one finds the 

general hydrodynamic equations: 

On . , . 

_ + V • (nv) = 0, (3a) 

( Q \ dP dll 

mn ^. +v . V y f + -^ +n — = , (3b ) 

de 

— + V • (ev) + V • Q + D^Pr = 0, (3c) 
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where we have defined the usual quantities: 

density: n(r, t) = J 7^l/( r > P> *)> ( 4 ) 

velocity: n(r,t)v(r,t) = J -^^^f(r,p,t), (5) 

pressure tensor : i%,( r > t) = m J -^0^ - ^ - v^j /(r, p, *), (6) 

energy density : e(r,t) = J ^ 2 ^ (p - mv) 2 /(r, p, t), (7) 

heat current : Q(r, t) = J ^ ~ ™v) 2 (j^ - v^j /(r, p, t), (8) 

1 / 3v Ov \ 

rate-of-strain tensor: D uv = - — — + — — - . (9) 

2 \dx v Ox,,; 

The first approximation for the distribution function is the local equilibrium form 

/ (0) (r,p,t) = {exp[/3(p - mv) 2 /2m + U - y\ - (10) 

where the thermodynamic variables /?, v and [i all depend on r and t, and U(r, t) has been 
defined after ([!]). This expression ensures that the collision integral (0) vanishes. This gives 
the lowest-order approximation inl/L, where I is the mean free path and L is a characteristic 
wavelength. If one uses ( |T0"D to calculate the quantities from (f|) to (H), one finds that Q = 
and 

n(r,t) = ±g 3/2 (z(r,t)), (11) 

P^(r, t) = 5^P(r, t), P(r, t) = ^g 5/2 (z(r, t)) = ~ £ (r, t). (12) 

Here z(r,t) = e ^t)M^)-u(r,t)] is the local f ugac it y , A( r ,t) = [2nh 2 /mk B T(r, t)] 1 / 2 is the 
local thermal de Broglie wavelength and g n (z) = J2iZi z l /l n are the well-known Bose-Einstein 
functions. The equilibrium value of the fugacity is given by zq{y) = e^ Q ^ l ~ u ^ : '\ with U{r) = 
Uo(r) + 2gno{v). Putting Q = and using (|T2"D in Eqs. (|3|), one obtains 

On 

— + V • (nv) = 0, (13a) 



mn 



<9v 



VP — nW, (13b) 



_ + | V (ev) = v ■ VP. (13c) 
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The linearization of the equations in ([13]) around equilibrium leads to the equations (4-6) in 
Ref. if we ignore the HF field in U(r,t). 

Solutions of these equations describe undamped oscillations, some example of which are 
discussed in Ref. || (see also Section IV). In order to obtain damping of the oscillations, we 
have to consider the deviation of the distribution function from the local equilibrium form 
(10). One assumes a solution of the quantum Boltzmann equation (1) of the form [|7,iT|: 



/(r, p, t) = /(°) (r, p, t) + / (0) (r, p, t) [1 + / W (r, p, f)ty(r, p, t), (14) 

where ip expresses a small deviation from local equilibrium. To first order in t/>, we can 
reduce the collision integral in (|2]) to 

2 f dp 2 f dp 3 f \x(p 2 ,P2 pI pI\ 

An9 J wrJ (^i rfP45(p+P2 - p3 - p4) H^ + ^"^"^J 

x/ (0) / 2 (0) (l + /i 0) )(l + /i 0) )(^3 + ^4 - ^2 - ll>) = L[lp], (15) 

where ipi = ^(r, p^, t). In the left hand side of ([!]), we approximate / by f^°\ The various 
derivatives of v(r, t),/j,(r, t),T(r, t) and U(r, t) with respect to r and t can be rewritten using 
the lowest-order hydrodynamic equations given in fll3|) . The resulting linearized equation 
for ip is (for details, see Appendix) 



u-VT 



2 



535/2 (; 



+ y^-^ ftt^ - l -8, v u 2 ) \ /<°>(1 + /W) = L[# (16) 



I T [2k B T 2g 3/2 (z)\ ' fe B T"^ 3' 
where the thermal velocity u is defined by mu = p — mv and the strain tensor is defined 
in (0). The linearized collision operator L is defined by (JTH). This equation can be shown 
to have a unique solution for ip if we impose the constraints 

J dpfM(l + / (0) )^ = / dpp^°\l + /<°>ty = / rfpp 2 / (0) (l + / (0) )^ = 0, (17) 

which mean physically that n, v and e [see (§),(§) and (|7|)] are determined only by the first 
term in (14). Since the left hand of ([16]) and also Z^ -* only depend on the relative 
thermal velocity u, the function if> will also depend only on u (i.e., not separately on p and 
v). 

It is convenient to introduce dimensionless velocity variables 

m \ x l 2 



2k B T 



u = t (18) 
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With these dimensionless velocity variables, ( |16l) becomes 

2 5 95/2{z 



71 

8a 2 mfc|T 2 



with = 
L'\ 



'2k B T\ 1/2 VT-£ 



m 



T 



2</3/2(«). 



+ 2Z^ ^ 



-<5 £ 2 



/( o) (0[1 + / (o )(0] 

= L'[V>], (19) 



l)- 1 and 



= / de a / ^3 / +Z2-t3- (e 2 + n - - ti) 

x/ (0) / 2 (0) (l + /f )(1 + /f)^ + ^4 - ^2 - ifi. 



(20) 



For a more detailed discussion of the mathematical structure of flT§| ) and (POD, we refer to 
the treatment of the analogous equations for classical gases (see for example, Ref. [12]). The 



most general solution of the integral equation (|Tj^) is of the form 



7T 



8a 2 mA;lT 2 

The functions and 5(0 obey the following integral equations: 

%5/2(2;) 



(21) 



l' mo] = € 



2^3/2(2;) 



/(o)(l + / (o)) 



^ - ^« 2 ) = (ex - l^e) / (o) (i + / (0) )- 



(22a) 
(22b) 



For (pi) to satisfy the constraints given in flTT|), we must also have 



/^ 2 ^(0/ (0) (i + / (0) ) = o. 



(23) 



Using a solution of the form (EH]) in conjunction with (|T^), one can calculate the heat current 
density Q in (§) and the pressure tensor in (g). One finds these have the form 



Q = -kVT, 

P»u = ~ 2^ 



D^v ~ g(TrD)«V 



(24a) 
(24b) 



where the last term in (|24b ) involves the non-equilibrium stress tensor. The thermal con- 
ductivity k and the shear viscosity coefficient ry are given in terms of the functions and 
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2kuT 



1/2 



48a 2 



m 



n 



m 
120a 2 



Hi 



1/2 



d& 4 A(O/ (0) (i + / (0) ), 

^4 jB(0/ (0) (l + / (0) ) _ 



(25a) 
(25b) 



Introducing (|24|) into the general hydrodynamic equations (§), one obtains the following 
hydrodynamic equations with the effect of shear viscosity and heat conduction included: 

dn 



dt 



+ V ■ (nv) 



0. 



d \ dP dU 
mn — + v • V + tj h n- — 



_d_ 

dx. 



277 



D^-i(T rJ D)o^ 



at 



V • (ev) + (V • v)P = V • (kVT) + 2r/ 



D 



fill 



(26a) 
(26b) 

(26c) 



We recall that n, v and e are still given by (4)-(7) with / = /(°), which means that the 
expressions in (11) and (12) for n,e and P are valid. The form of these equations can be 
shown to agree with those originally obtained by Uehling and Uhlenbeck J7|. We note that r\ 
and k are position-dependent, but only through their dependence on the equilibrium value 
of the fugacity z = Zo(r). One slight generalization we have made over the derivation in 
Ref. |7j is that we have included the Hartree-Fock mean field. 



III. THE TRANSPORT COEFFICIENTS 
A. The thermal conductivity 

To find the thermal conductivity as given by (|25a|) , we can introduce a simple ansatz for 
the form of the function A(£) g[T2|: 



2 %5/2^) 



(27) 



2^3/2(2) 

The constant A is determined by multiplying ( |22a| ) by £[£ 2 — 5g 5 /2(z)/2g 3 / 2 (z)] and inte- 
grating over 



A— J d^e 



595/2(2) 
2^3/2(2) 



+ /(°); 



505/2 (*)' 
203/2W, 



5^5/2 jz) 
203/2(2;) 
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15vnV2 
4/a 



7 , , 5^ /2 (z) 
2*»« - 2^) 



(28) 



where the integral is defined by 



^2 _ %5/2(» 



2 _ 5fl5/2(*) \ ' 



2^3/2 (^) 



(29) 



In order to evaluate the integral in (p9"D , it is convenient to introduce the change of variables 



(30) 



Then we introduce transformations from (,' x £.yC z to and CxCyCz to £ where 9', 9" 

and 0', 0" are the polar and azimuthal angles with respect to the vector £ . One obtains the 
following expression for I a-, 



I A = -4v / 2tt 3 /' 



Ai 



f>o f>' Cdyi f 1 dtftifFfa, e,y',y";z)[y" + y" 2 -2y' 2 y"X (31) 
Jo Jo Jo Jo 



where y' = cos 9', y" = cos 9". Here the function F is defined by 



F = / (o) / (0) (1 + /f)(1 + / (o )) 



(1 - ze-?)(l - ze~^)(l - ze -«)(l - ze-#) 



(32) 



with 



e 3 = \(e + 2&ry" + e' 2 ), e 4 2 = ^ 2 - + r 2 )- 



(33) 



Inserting the expression in (|27D into fl25a|) and carrying out the integration, we obtain 
the following expression for the thermal conductivity k: 



75kn fmknT\ 1/2 vr 1 / 2 



Q4a 2 m \ tt 



iQl'Az) 



7 , , 59 2 / 2 ( z ) 



(34) 



where the function I' A {z) is defined in (|31"|). 
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B. The shear viscosity 



In evaluating the shear viscosity in (|25b|) , the simplest consistent approximation |f||l2| is 
to use -£>(£) = B. The constant B can be determined by multiplying ( |22b| ) by (C^— ^C 2 /3) 
and integrating over £, 



6 = ^ /^(^-^« 2 ) 2 / (0) (i + / (0) ; 



i 

3' 



5^5/2 ( 



2 



2/ 



a; 



The function I B is defined by 



Ib= I dZ[^ u -~-5, v e)L' 



1 

— < 

3 



= -2v / 27r 3 / / f 



(35) 



(36) 



where 

POO POO pi pi _ 

is W = / / # / dy> / dy"F(e , y', y"; *)£o 2 £ /7 (i + y' 2 + y" 2 - Sy'V 2 )- (37) 

Jo Jo Jo Jo 
This involves the same function F as defined in (132]) . Using the ansatz = -B in 

conjunction with (|35|) in ( 25b ), we can carry out the integral. Our final the expression for 
the viscosity 77 is 

5 1 (mk B T\ l/2 ir 1 / 2 



16 a z 



TV 



W) 95/2{z) ' 



(38) 



C. High-temperature limit 

The formulas in (|34]) and (^) give « and rj as a function of the fugacity 2. We note 
that the integrals V A and I' B also depend on z. These four- dimensional integrals in (|31|) and 
(pTD can be evaluated numerically. In Figs. 1 and 2, we plot both k and 77 as a function of 
the fugacity 2. These graphs are valid for any trapping potential Uo(r) since the latter only 
enters into the equilibrium fugacity z = exp[/3o(/x — ^o( r ) — ^9 n o( T ))]- 

In the high-temperature limit, where the fugacity z = e^- u W is small, the local dis- 
tribution function can be expanded in terms of z. To third order in z, the function F 
in (]3~2D reduces to 
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e -(ti+n + 2z 3 e-^ + ^[cosh(^'y') + cosh^Y')] + 0(z 4 ). 



(39) 



The integrals in (^T|) and 

r B (z) = 2i' A ( z ) 



l + z- 



(40) 



can be evaluated analytically, 

9 /3 
16 V 2 

We thus obtain the following explicit expressions for the transport coefficients to first order 
in the fugacity z: 

7^2 9 



1 /75\ k-Q ( k^T 

K = — 



8 V64/ a s 

1 (75\ kB (knT 
8 V647 a 2 



1 / 5 \ m ( knT 



V 8 U6/ a 2 

_ 1 / 5 \ m fk B T 
~ 8 V167 a 2 




(41a) 



(41b) 



The local equilibrium density n (r) in the high-temperature limit is the classical result 
tiq = A~ 3 e^^~ u ^ . The terms first order in z in (41a) and (|41b| ) give the first order 
corrections to the classical results due to Bose statistics. If we ignore the HF mean field 
2gn (r), these transport coefficients reduce to the expressions first obtained for a uniform 
Bose gas by Uehling ||. The Bose quantum corrections to the classical results in both r\ 
and k depend on the local fugacity z and, through this, on position in a trapped gas. 

The density-independent terms in Eqs. (|41a ) and ( |41b|) are 8 times smaller than the well- 
known Chapman-Enskog expressions for classical hard spheres. This is due to the difference 
(see also Ref. [H) in the quantum binary scattering cross-section for Bosons when correctly 
calculated using symmetrized wavefunctions (a = 8iia 2 , instead of na 2 ). 



IV. HYDRODYNAMIC DAMPING OF NORMAL MODES 



The linearized version of the hydrodynamic equations in (p6|) are 
d5n 



dt 



+ V • (n 5v) = 0, 



mriQ- 
85P 



d5v. 



d5P . dU n d5n d 
— on- ZgnQ— h 



dx. 



dx. 



dx n dx u 



27] 



D/iu - -(TrD)^ 



Jv • (P 8v) + j^v ■ VP + |v • (kVT), 



(42a) 
(42b) 

(42c) 
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where repeated Greek subscripts are summed over. 

Taking the time derivative of ( [42 b| ) and using ( [42 a| ) and ( |42c| ) gives an equation for the 
velocity fluctuations which is only coupled to the temperature fluctuations: 

m d 2 5v, = 5Po _ _d_ _ 2 

<9t 2 3n <9x M ftr^ 3 <9a^ 

+ 2 ^S(V • Sv)] - ^(V • *v)^-no 



1 5 In 



dt ^ 3 I dt ^ 



lk^ V -^ 6T) - (43) 



Here we have used the relation VPo = —no'VU valid in equilibrium. Since we assume the 
dissipative terms to be small, we can use the lowest-order hydrodynamic equation for ST in 
(H) (see (A7) in Appendix), 

^ = - 2 - Tq V ■ Sv. (44) 
dt 3 v ; 

Equations ( fl3"D and (0) are thus a closed set of equations for the fluctuations. Equation fl4"3"|) 
generalizes that given by (13) in Ref. to include damping due to viscosity and thermal 
conductivity as well as the HF mean field 2gn. 

At this point, we introduce a crucial aspect of a trapped gas which leads to an important 
modification of the preceding hydrodynamic equations. Kavoulakis et al. |J have pointed out 
the breakdown of the local equilibrium (hydrodynamic) description in the low density outer 
region of a trapped gas. The corrections involving thermal conduction and viscous processes 
are expected to vanish (over a mean free path) in the dilute region where collisions become 
ineffective in producing local equilibrium. This effect can be simulated in the hydrodynamic 
equations by multiplying the hydrodynamic expression for the thermal conductivity in ( |2 la] ) 
and the viscous stress tensor in ( |24b| ) by a step function that vanishes outside the cloud : 



Q = _ K VT0(r o -r), 



Dfj,u — -(TrD)Sfj, u 



G(r - r), 



(45a) 
(45b) 



For simplicity of notation, we only consider an isotropic trap, in which case the hydrodynamic 
description breaks down on a sphere of radius r . A method of calculating r is given 
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in Ref. || for a classical trapped gas. The equation of motion in ( fl"3|) for the velocity 
fluctuations is modified to 



d 2 5v, 



5P n d 



ni- 



di 2 



(V ■ <Jv) 



dx 



_( 5v . V C/ )--(V-5v)^ 



3n 9^ 

+2^ [no(V . ( 5v)]-^(V-5v) 9//i: 



dx h 

1 9 f 
+ — — <2 V 



n dx v 
2 



9* Mly 3 I dt ^ 



6(r - r) 



V- [KV8TQ(r -r)]. 



(46) 



3n 9x M 

We now look for normal mode solutions of (|46|) with the convention 5v(r, £) 
5v w (r)e-^*, ST(r,t) = <JT w (r)e _<w * and £>(r,t) = D^e - *"*. Combining (||) and (|6 
we obtain a closed equation for ^v^: 



-mn (u 2 8v, 



d 



—iuj- 



2r] 



9(r - r) 



9n u dx^ 

where the operator W^v^] is defined by 



dx v 

^^^-{V • [kV(V • 5vJ6(r - r)]}, 



(47) 



5P d d 2 ,0£/ o 

V • ov,„) - - — (dv u ■ VU ) - ~(V • dv u 



dx. 



3n dx^ 

+2g^-[n (V ■ 6v u )) - i?(V • Sv^"" 



<9x, 



(48) 



The undamped solution dv^ of ( f47j ) when k, 77 = is given by Wjt[5v w ] = Q 2 Sv u! ^ (a tilde 
denotes undamped solution). For later use, we note that the operator has the following 
property 

(49) 



noSvu^W^Sv^dr = / n 5v LUI1 W ll [5v^}dr. 



Therefore the undamped solutions satisfy the following orthogonality relation 



noSvl, W^Sw^dr = 0, if Q' ^ Co. 



(50) 



We consider an isotropic trap Uo(r) = ^vnio^j . The undamped solution for a monopole 
or breathing mode is described by 5v w (r) ~ r, the undamped frequency is found to be 



12 



uj = 2lu [[|. Using this solution in (|47j), we find D uj^i/ ~ & jiv and hence there is no shear 
viscous contribution to the damping of this mode. In addition, since V • <5v^ =const, the 
last term in ([IT]) involving the thermal conductivity coefficient k does not contribute. We 
conclude the monopole mode has no damping for an isotropic harmonic trap. As noted in 
Ref. 0, this surprising result was first derived for a classical gas by Boltzmann in 1876. 

The undamped solution for divergence-free surface modes (V • <5v w = 0) is given by |3| 
<5vu,(r) = Vxa,(r), where 



X w (r) oc r l Y lm (9,(j)), 



(51) 



with the dispersion relation u = y/loooi I — 1, 2, ■ ■ -. Using the solution (|51~D in the right 
hand side of (0), we find (there is no contribution from the thermal conductivity because 
V-5v w = 0). 



RHS of (47) = -2ito 



_d_ 

dx,. 



9(r -r) 



2iu 



z dri 2 77 . 

— — — muj e(r -r) + -6(r- r ) 
KbT oz r 



2iu(l - 1) 



z dri 2 . , ri , . 
— — — mw o e(r - r) H b(r - r ) 

/vB-t OZ T o 



5v, 



(52) 



For I = 1, we see that the right hand side of (47) vanishes. This means that the center- 
of-mass mode with frequency Uq has no hydrodynamic damping. For I > 1, the undamped 
solution is no longer the solution of the equation of motion In this case, Eq. ( f4T|) 

must be solved for r < r with a boundary condition at r = r to take into account the 
discontinuity at r given in (p2j). 

To complete this section, we derive a general expression for hydrodynamic damping 
from the linearized hydrodynamic equations for the velocity fluctuations given by (fj7p. 
Multiplying (|47D by <5t> * M and integrating over r, we obtain 



-mio 



n \Sv u \ 2 dr + m n 6v* W^Sv^dr 



2i] 



D 



LUfll/ 



0(r o — r) j dx 



4T 
9noo; 



2r/ 



UJ/1U 



-(TrDJV 



4T 

+ ?r ^|V(V-5v a 

9nour 



6(r — r)cfr. 



(53) 
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Since the undamped solutions ^v^ satisfy the orthogonality relation (50), one sees that these 
undamped solutions can be used to find the eigenvalue uj 2 of ( |53D to first order in r\ and k. 
Inserting the undamped solutions into fl5*B|), we find it reduces 



m{pj — uj ) no\8vu\ dr 



2q 



+ 



4T n 



6(r — r)dr. 



(54) 



9n uj 2 

Assuming the complex solution uj = uj — iT, one finds the damping rate T is given by the 
expression (valid to first order in rj and k) 



2r] 



4T 
9n u 2 



0(ro — r)dr 



2m I noldv^dr 



(55) 



The damping rate given by (|55f) is precisely 0] the one obtained in Ref. 0, where the 



damping rate is calculated from the rate of increase of the total entropy [14]. The transport 
coefficients k and 7] in (|55| ) are given in Figs. 1 and 2 for a trapped Bose gas and in general 
are dependent on position through the fugacity. 



V. GENERALIZATION TO A SUPERFLUID GAS 

In this section, we extend the hydrodynamic equations derived in Section II to the case 
of a superfluid trapped gas by generalizing the results of Ref. [|J to include hydrodynamic 
damping. The two-fluid equations in Ref. consist of equations of motion for the conden- 
sate and a set of hydrodynamic equations for the non-condensate. The latter description 
is only valid for in the semiclassical region k^T ^> huo,gno. The time-dependent Hartree- 
Fock- Popov equation of motion for the condensate wavefunction can be rewritten in terms of 
a pair of hydrodynamic equations for the condensate density n c and the superfluid velocity 



771 



dn c 
dt 



-V ■ (n c v s ), 



(56a) 
(56b) 



The potential 4> is defined by 
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h 2 V 2 \n (r tY\ l l 2 

<Kt, t) - ~ 2m[n c (r,t)] V* + + ^ f) + ^ ^ (5?) 

where n(r, £) is the non-condensate density describing the excited atoms. The first term in 



( p7|) is sometimes called the quantum pressure term and vanishes in uniform systems. It is 
associated with the spatially varying amplitude of the Bose order parameter. 

The equations of motion for the non-condensate can be derived from the kinetic equation 
(H) for the distribution function /(r, p, t) of the excited atoms. The HF mean field in 
U(r,t) is now 2gn(r,t) = 2g[n(r, t) + n c (r,t)], where n(n c ) represents the local density of 
the non-condensate (condensate) atoms. Apart from this change, the expression for the 
collision integral in (0) is still applicable at finite temperatures below Tbec- Thus the 
hydrodynamic equations for the non-condensate can be derived using precisely the same 
procedure developed in Section II for T > T B ec- We obtain equations analogous to ( PB|) in 
the following form: 

OTl 

— + V • (nv n ) = 0, (58a) 
at 



d _\ dP dU d 



m ,^_ + v „.Vj tV , _ 

de 

— + V ■ (ev„) + (V ■ v n )P = V • («VT) + 2 V 



du d r r i 



- ~(TrD)5„ v 



(58b) 
(58c) 



Here n, P and e are given by the same expressions as n, P and e in (^) and (0). D^ v is 
defined as in @, where now the velocity v is the normal fluid velocity v n . The expressions 
for the transport coefficients given in (j34D and (|38|) are directly applicable below T B ec f° r 
any trapping potential, apart from change in the form of the HF field in the local fugacity 
z. As in the case above Tbec considered in Section IV, one must introduce a cutoff in the 
linearized hydrodynamic equations to take into account that the terms proportional to the 
transport coefficients vanish when the density becomes too low to sustain local equilibrium. 



We might note that the results for k and r\ by Kirkpatrick and Dorfman [11,15] for a 
uniform gas below Tbec have an additional factor (1 + 6n c A 3 ) -1 , where b is of order 1. This 
factor arises from collisions between condensate and excited atoms, which are important at 
very low temperatures (T <C Tbec)- Our present analysis ignores this contribution. In this 
regard, we note that (|1|) and (|2|) are based on validity of a simple semi-classical Hartree-Fock 
description of excited atoms (ignoring the off-diagonal or anomalous self-energy terms). This 
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has been shown [|16j to give an excellent approximation for the thermodynamic properties 
of a trapped Bose gas down to much lower temperatures (T <C T BE c) than in a uniform gas. 
If we neglect the first term in (|57|) , the equations in (|5^| ) and can be used to derive 



JD the two-fluid equations for a trapped Bose-condensed gas in the Landau form JET], except 
that our equations do not include any effect from the bulk viscosities. In a uniform superfluid 
Bose gas, Kirkpatrick and Dorfman [IT] found that the bulk viscosities vanish in the finite 



temperature region (k^T ^> grio) we have been considering. Our present results for a trapped 
Bose gas are consistent with this and moreover justify the omission of any dissipative terms 
in the condensate equation of motion (|56b| ). 

VI. CONCLUDING REMARKS 

Summarizing our main results, starting from the microscopic equations of motion, we 
have derived a closed equation (|43|) for the velocity fluctuations of a trapped Bose gas, 
which includes damping due to hydrodynamic processes. We have also obtained explicit 
expressions for the shear viscosity r\ and the thermal conductivity k for a trapped Bose gas, 
which depend on position through the local fugacity z . We have given a detailed derivation 
for a trapped gas above Tbec, and more briefly discussed the generalization of these results 
to a Bose-condensed gas below Tbec in Section V. 

For illustration, in Section IV, we used the linearized hydrodynamic equations to discuss 
some of the normal modes of oscillation for a normal Bose gas above Tbec- We show that 
the monopole mode of an isotropic trap was undamped. We also note that our linearized 
equations in ( |42b[ ) and ( f42c| ) do not take into account that in the low density region in 
trapped Bose gases, the hydrodynamic description breaks down, as emphasized in Refs. ||[| . 
We show that if a spatial cutoff is used in our linearized hydrodynamic equations, the 
damping of the normal mode solutions is given by the same expression as one finds from 
calculating the rate of entropy production |||||. In this formalism, one has reduced the 
calculation of the damping to that of determining the boundary where hydrodynamics breaks 
down. This has been worked out in Ref. for a classical trapped gas, in which the transport 
coefficients are independent of density. An interesting problem for the future is to use the 
semi-classical kinetic equation given by (H) and (0) to give a more microscopic analysis of 
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the low density region where expanding around the local equilibrium solution [see (14)] is 
no longer valid. 

We note that the general form of the hydrodynamic equations given in Eq. (PB| ) can 
be derived more generally using conservation laws for conserved quantities in conjunction 



with local thermal equilibrium. This approach is developed in many textbooks [18| and 



can be formally extended to deal with inhomogeneous superfTuids in an external potential 



(see, for example, Ref. ||19|| ). This method is, of course, phenomeno logical in that it in- 
troduces various linear response coefficients [relating densities (thermodynamic derivatives) 
and currents (transport coefficients) to fields] whose evaluation requires some specific mi- 
croscopic model. While very useful, this approach must be ultimately justified by a fully 
microscopic treatment given by the underlying kinetic equations, such as used in this pa- 



per and in Refs. [0,||,0,|TTJ. In particular, as noted above, such a microscopic approach is 



needed to deal with the breakdown in the local equilibrium solution in the low density tail 
of trapped gases. Moreover we note that our microscopic derivation in Section V (see also 
Ref. |Q) gives a more complete description, as can be seen from the fact that we obtain 
separate conservation laws for the condensate and non-condensate densities. The standard 
phenomenological two-fluid equations do not determine the condensate and non-condensate 
density fluctuations separately, but only the total density fluctuation. 
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APPENDIX A: 



We briefly sketch the derivation of the left hand side of the kinetic equation in ( p"6| ) 
Using ( |T0| ) in the left hand side of ([I]), one has 
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d p 

dt m 



V r U{r,t)-V p 

,2 



/ (0) (r,P,t) 



1 f d p \ m« 
z 1 eft m r I 2k B T 2 \dt'm 



4 + ^V)t 



mu / 9 p _ \ V r U(r, t) 

' + — • V r v + r , „ 7 ■ u 



/(«)(! + /W). 



(0)N 



(Al) 



/cb^ \dt m J k B T 

Using the lowest-order hydro dynamic equations in flT3"|) and the expressions for the density 
n in (|TTD and the pressure P in (|12D , one finds 



dn 3n dT jk-sT dz 
~dt ~ 2T~dt + z dt 



3n m 

= — -v • vr 

2T 



lk B T 



v • Vz — nV • v, 



where 7 = ^ jsgi/ 2 (z), and 



dP _5PdT nk B T dz 
~dt ~ 2T~dt + z dt 

= — — v • Vi v ■ Vz P V • v 

2T z 3 v ' 



One may combine these equations to obtain 

-v • Vz = 



<9z 



-i-v-fvP-^VT 

nfc B T V 2 T 



(A2) 
(A3) 



(A4) 
(A5) 



(A6) 



dT 
~dt 



-TV ■ v - (v • V)T. 



(A7) 



The analogous equation for dv/dt is given directly by (13b). Using these results in (Al), 
one finds that it reduces to 



9 p _ 

dt m 



V r U ■ V p j f 

,2 5p 



(0) 



iu-VTf^l 

T V 2fc BT 2nfc B T 



u 



u • (u • V)v V ■ v 



/<°>(l + / (0) ), (A8) 



where we recall u = p/m — v. This can be rewritten in the form shown on the left hand 
side of (0). 
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FIGURE CAPTIONS 



Fig.l: The thermal conductivity k as a function of the fugacity z. The values are normalized 
to the classical gas results at z = 0. (see ( |41a]j ) 

Fig.2: The viscosity coefficient rj as a function of the fugacity z. The values are normalized 
to the classical gas results at z — 0. (see ( [41 bp ) 
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Figure 1 
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